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We study the static screening in a Hubbard-like model using quantum Monte Carlo. We find 
that the random phase approximation is surprisingly accurate almost up to the Mott transition. We 
argue that in alkali-doped Fullerenes the Coulomb pseudopotential fi* is not very much reduced by 
retardation efTects. Therefore efficient screening is important in reducing fi* sufficiently to allow for 
an electron-phonon driven superconductivity. In this way the FuUerides differ from the conventional 
picture, where retardation effects play a major role in reducing the electron-electron repulsion. 
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The random phase approximation (RPA) has been 
very widely used in solid state physics. It properly de- 
scribes the screening when the kinetic energy is much 
larger than the interaction energy.EI In the opposite limit, 
however, the RPA is qualitatively wrong. Little is known 
about the more interesting situation when the two en- 
ergies are comparable. In this paper we show for a 
Hubbard-like model that the RPA gives a surprisingly 
accurate description of the static screening on the metal- 
lic side of a Mott transition until the system is close to 
the transition. 

For conventional superconductors the electron-phonon 
interaction leads to an effective electron-electron attrac- 
tion. This interaction is counteracted by the strong 
Coulomb repulsion, which is, however- Jpelieved to be 
strongly reduced by retardation effects. 0'cl The resulting 
effective Coulomb interaction is described by the dimen- 
sionless Coulomb pseudopotential /i*, which is believed 
to be typically of the order 0.1. Here we argue that the 
situation for A3C60 (A= K, Rb) is different. We find 
that retardation effects are rather inefficient. Therefore 
the screening of the Coulomb interaction becomes impor- 
tant for reducing the electron-electron repulsion. Thus, 
although the superconductivity|-|in A3C60 is driven by 
the electron-phonon intcraction,Q the origin of the strong 
reduction of /i* is different from the current picture of 
conventional superconductors. In the scenario we are 
putting forward, several puzzling phenomena find a nat- 
ural explanation. In A3C60 (A= K, ^b) the transition 
temperature Tc is reduced by pressure^ For CssCgo, how- 
ever, which only under pressurci-becomes a superconduc- 
tor, Tc increases with pressure.Q This is consistent with 
the picture where /i* is reduced by screening, since the 
screening is less efficient close to a Mott transition. Sec- 
ond, it was very early pointed out that the alkaJi phonons 
ought to couple efficiently to the electronsjj althoug 
later experiments showed that this was not the case 
We show that efficient screening reduces the coupling to 
the alkali phonons. 

We first discuss the screening in the RPA. In the ran- 
dom phase approximation it only costs kinetic energy to 
screen a test charge. In the limit where a typical Coulomb 



integral U is large compared with the band width W, the 
kinetic energy cost of screening is relatively small com- 
pared with the potential energy gain, so the screening is 
efficient. This means that as a test charge q is introduced 
on a site c, almost the same amount of electronic charge 
moves away from the site, leaving it almost neutral. This 
argument neglects, however, that when an electron leaves 
a site it has to find another site with a missing electron or 
there is a large Coulomb energy penalty. Thus the RPA 
is accurate for small values of U/W, while it is qualita- 
tively wrong for large values. It is not clear what happens 
for intermediate values. 

To study the screening in A3C60, we use a Hubbard-like 
model, including the three-fold degenerate ti„ orbital: 

^= E E hm.jm'^l^, ^jm'a ^^) 

+ UJ2 J2 n^raan^m^a'+<iU^ncma■ 

i (mcr)<(m'(T') "^^ 

The first term describes the kinetic energy, the second 
term the on-site Coulomb interaction and the third term 
the interaction with the test charge q on site c. ipim,cr 
annihilates an electron on site i with orbital quantum 
number m and spin a, andj|U^cr = V'LaV'jmo-- The ef- 
fect of orientational disordeiolEl is built into the hopping 
integrals t^^^^' . The band width is about 0.63 eV. Mul- 
tiplct effects are not included, but we remark that they 
tend to be counteracted by the Jahn- Teller effect which 
is also neglected. The test charge is assumed to interact 
with the electrons on the same site via the Coulomb in- 
tegral U. The system has three electrons per molecule, 
i.e. a half-filled ii„ band. 

We have investigated the model by using a_La±±ice dif- 
fusion quantum Monte Carlo (QMC) method. EiJll3 In this 
method a trial function \4't) is constructed and allowed 
to diffuse towards the exact solution, under the constraint 
of a fixed node approximatioi!.-, \iPt) is obtained from a 
generalized Gutzwillcr Ansatztil 



I^t) 



D ric 

9 5o 



*o), 



(2) 



where |^o) is a Slater determinant constructed from so- 
lutions of eqn. (|l|) in the Hartree approximation, D is 
the number of double occupancies in the system, and ric 
is the number of electrons on site c. g and go are vari- 
ational parameters, g^ is the usual Gutzwiller factor 
while gg" allows us to optimize the charge on the site 
with the test charge. In addition to the DMC calcula- 
tion, we also perform a variational Monte Carlo calcula- 
tion (VMCI_and the energy is minimized as a function 
of g and go S3 In all cases the state is assumed to be para- 
magnetic. For U/W ~ 2.5 ther&is a transition to an an- 
tiferromagnetic Mott insulator)13 where the screening is 
very inefficient. Here, however, we focus on U/W < 2.5. 
We obtain the charge on site c from the extrapolated es- 
timator He « 2nc{DMC) - ndVMC), where ndVMC) 
is the expectation value for the wave function (0) cal- 
culated by VMC and nc{DMC) is the mixed estimator 
from the DMC calculation. 

To test the accuracy of the approach, which involves 
the fixed-node approximation and uses the extrapolated 
estimator, we have compared the results of our QMC 
calculations with the the exact results from exact diago- 
nalization of a system with four molecules (12 electrons). 
The comparison shown in Fig. |l| illustrates that the QMC 
calculations are quite accurate for the system we are an- 
alyzing here. 

Since we are interested in the linear response, we 
should calculate the effect of an infinitesimally small test 
charge q. Because of the statistical error in a QMC calcu- 
lation it is, however, difficult to determine the response 
to a small perturbation. To get a good signal-to-noise ra- 
tio, we would therefore like to use as large a test charge 
as possible. To estimate how large we can make q and 
still be in the linear response regime, we have performed 
Lanczos calculations for a range of different test charges. 
We find that for q < 0.25 e the response is practically 
linear. 
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FIG. 1. Screening charge An on the site of the test charge 
{q — 0.25 e) as a function of U/W, where U is the Coulomb 
interaction and W is the band width. Exact diagonalization 
and QMC calculations have been performed for four molecules 
(12 electrons). The figure shows that the QMC calculations 
are quite accurate over the whole range of U/W . 
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FIG. 2. Screening charge An on the site of the test charge 
(q = 0.25 e) as a function of U/W, extrapolated to infinite 
cluster size. The full curve shows the screening charge in the 
RPA, obtained from Hartree calculations for the Hamiltonian 
(hi). The crosses with errorbars give the results of the QMC 
calculations. The RPA screening remains rather accurate up 
to U/W ~ 2, but fails badly for larger values of U/W. The 
screening is very efficient for U/W ~ 0.5 — 2.0. 

We have performed QMC calculations for larger clus- 
ters of A^moi= 32, 48, 64, 72, and 108 molecules, where 
exact diagonalization is not possible. The screening 
charge An^ — nc{0) — ndq) was extrapolated to infi- 
nite cluster size, assuming a finite-size scaling of the form 
Anc(Nmo\) — Arij, -I- a/Nj^^i- The results are shown in 
Fig. |. For rather small values of U/W (- 0.5 - 1.0), 
the RPA somewhat underestimates the screeruHg. Such 
a behavior is also found in the electron gasJl3 For in- 
termediate values of U/W (-- 1.0 - 2.0) the RPA gives 
surprisingly accurate results. This is one of the main 
results of this paper. For large U/W, the RPA rapidly 
becomes qualitatively wrong, as discussed earlier. We are 
now in the position of addressing the superconductivity 
in AsCeo- 

In the theory of superconductivity, a dimensionless 
quantity fi* , the Coulomb pseudopotential, is introduced 
to describe the effects of the Coulomb repulsion. One 
introduces /i = UN{0), where U is a typical screened 
Coulomb interaction and N{0) is the density of states 
per spin at the Fermi energy. Retardation effects renor- 
malize fi to fj,*, and are described by laddetxiiagrams in 
the statically screened Coulomb interaction.aij It is found 
that 
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where Ugi and tUph are typical electron and phonon energy 
scales, respectively. Often fAn(uJei/i^ph) is substantially 
larger than unity. In that limit the last part of eqn. (0) 
holds, i.e. fi* is determined solely by retardation effects, 
independently of the screening, which only changes fi. 

In solid Cgo we have many narrow subbands [^ 0.5 eV 



wide) , spread over a range of about 30 eV. In the tra- 
ditional approach one assumes that the relevant energy 
range extends over all this region. Summing the ladder 
diagrams in the screened Coulomb interaction leads to a 
large renormalization of fi. Exact results for a two-band 
model show, however, that in the appropriate limits this 
approach greatly overeaidmates the renormalization due 
to the upper sub band£3 In the limit when a sub band 
is far away from the Fermi energy, the correct approach 
is to first project out the high energy degrees of freedom 
corresponding to this sub band. This leads to an effec- 
tive Hamiltonian, expressed in terms of the unscreened 
Coulomb matrix elements, which describes the low en- 
ergy properties of the system. The main difference be- 
tween the two approaches is the order in which high and 
low energy degrees of freedom are treated. In the tra- 
ditional approach the Coulomb interaction is screened 
first, which in particular involves the low energy degrees 
of freedom. After this the high energy degrees of freedom 
are projected out. This approach involves uncontrolled 
approximations. Our approach, instead, projects out the 
high energy degrees of freedom first, and it allows us to 
make statements about the importance of these degrees 
of freedom. Although these arguments were presented in 
the context of Ceo, they are rather general. We now make 
more specific arguments for Ceo to provide further evi- 
dence that the retardation effects from higher sub bands 
are not very large. 

From Auger measurements on KeCeo, the Coulomb in- 
teraction U between two holes in an otherwise full ti„ 
band has been estimated ta, about 1.5 eVO This reduc- 
tion of U from about 4 eVEfl for a free molecule to about 
t^insui = 1-5 eV for the the insulating solid, is mainly 
due to intramolecular processes and to polarization of the 
molecules surrounding the two holes. Since the excitation 
energy of the relevant final state in the Auger experiment 
is rather small (about 1.5 eV), C/insui should contain the 
renormalization from all the higher sub bands, except 
possibly the ones closest to the ii„ band. If we multiply 
^insui by A^(0) ^ 6,E3 the result is a very large /^ ^ 9, 
much too large to allow for a phonon induced supercon- 
ductivity unless fi is further reduced by other effects. 

In KaCeo screening and retardation effects inside the 
tiu band become available. The argument against sum- 
ming ladder diagrams in the screened interaction were 
only justified for higher sub bands. Within tfagiiiu band 
we therefore rely on this conventional theory,EHj which in 
addition usually uses Thomas- Fermi or RPA screening. A 
priori, the use of RPA seems highly questionable for these 
strongly correlated systems. Our calculations, however, 
support this approximation! unless the system is close 
to a Mott transition. Taking the long range Coulomb 
interaction inio account, the RPA screening reduces /i 
to about 0.4.I1j Including the additional retardation ef- 
fects inside the tiu band according to eqn. dsh finally 
renormalizcs /i to /i* « 0.3. Thus the Coulomb pseu- 
dopotential is primarily reduced by screening and not 
by retardation effects. In contrast, using eqn. (W with 



ujei ~ 15 eV and ujph ~ 0.1 eV would result in fi* < 0.2, 
practically independent of fj,. A Coulomb pseudopoten- 
tial /i* « 0.3 is substantially larger than for conventional 
superconductors,D but it is not so large that it prevents 
the superconductiv^ from being driven by the electron- 
phonon interaction.Ej Recent tunneling experiments give 
/I* = 0.329 for RbgCeo.B 

We now turn to the question how Tc changes with the 
lattice constant a. The main effect of increasing a is to 
decrease the band-width W and increase the density of 
states at the Fermi level N{0). Using McMillan's for- 
mula, Tc is given by 



Tr 



T2 



exp 



-1.04(1 + A) 



A-^*(1 + 0.62A) 



(4) 



with A = N{Q)V the electron-phonon coupling constant. 
^* is calculated from ^ — N{Q)U\nsu\ (1 — 7), where C/insui 
is a typical unscreened Coulomb matrix element and 
7 = dn/dq describes the screening within the ti„ band. 
Assuming that ujei in Eq- @ is large, /j,* is practically in- 
dependent of the lattice constant a. Since iV(0) increases 
with decreasing a, the electrattphonon coupling A be- 
comes stronger, increasing TcE3 Assuming a small We/, 
corresponding to the ii„ band width, it is no longer true 
that /i* is independent of /x. However, if the RPA is valid, 
/i is almost independent of the lattice constant, since the 
increase in iV(0) is counteracted by a slightly more ef- 
ficient screening 7i?pA (cf. Fig. 0). Hence also in this 
scenario we find that Tc increases with a. But what hap- 
pens when the lattice constant a becomes large enough 
that we enter the region where the screening starts to 
break down? Then /i will start to increase considerably 
with a. Assuming a large ujd, ^* is still independent of 
/i, and therefore Tc should keep increasing. For small 
ujci, on the other hand, /i* will start to rapidly increase 
with a, leading to a steep drop in Tc- This resembles 
the anomalous behavior observed in CsaCeo: it only be- 
comes superconducting under pressure, with Tc rapidly 
decreasing with increasing lattice constant .El 

It might appear that efficient screening is not re- 
ally helpful for superconductivity. Phonons couple to 
the electKjMS by perturbing the potential seen by the 
electrons.0i3 An example being the longitudinal modes 
of a jellium. Efficient screening tends to weaken the 
coupling to such phonons, since it reduces the pertur- 
bation considerably. To some extent, such a reduction 
also seems to be at work in Ceo. Initially it was ex- 
pected thatjihe coupling to the alkali phonons would be 
very strong.!] Each Ceo molecule is surrounded by 14 al- 
kali ions with relatively weak force constants. When an 
electron arrives on a Ceo molecule one would therefore 
expect that the surrounding alkali ions respond strongly. 
This was, however, not confirmed by experiment. For 
instance, an alkali isotope effect wuld not be observed 
within the experimental accuracy.Q This finding can be 
naturally understood as an effect of the efficient screen- 
ing found in our calculations. When an electron arrives 



on a Ceo molecule, other electrons leave the molecule, 
which thus stays almost neutral. The alkali ions then 
only see a small change in the net charge and therefore 
couple weakly. In a similar way it follows that intramolec- 
ular phonons of Ag symmetry couple weakly. An Ag 
phonon shifts all the tiu levels on a given molecule in the 
same direction. This shift of the center of gravity can be 
screened very efficiently by transferring charge from the 
molecules where the levels move upwards to those where 
they move downwards. The modes that are important 
for the superconductivity in solid Cgo are, however, dif- 
ferent. An intramolecular Hg phonon does not shift the 
center of gravity of the ii„ level. Thus the Hg phonons 
are not screened by the transfer of charge. Hence for 
these phonons the efficient screening serves to reduce /i* 
without affecting the electron-phonon coupling. 

To summarize, we have calculated the static screening 
of a point charge for a Hubbard-like model using quan- 
tum Monte Carlo. We find that the RPA is surprisingly 
accurate up to values of U/W fairly close to the Mott 
transition. For larger U/W the screening rapidly breaks 
down. This result should have quite general implications 
for the physics of systems close to a Mott transition. Here 
we have studied the consequences for the superconduc- 
tivity in the alkali-doped FuUerenes. We have provided 
arguments that for AaCgo (A= K, Rb) retardation ef- 
fects are very inefficient in reducing the electron-electron 
repulsion. Instead, and unlike for textbook supercon- 
ductors, screening is mainly responsible for the reduc- 
tion of the Coulomb pseudopotential fi* . This results in 
a n* small enough that the electron-phonon interaction 
can drive the superconductivity. Nevertheless /i* is sub- 
stantially larger than for conventional superconductors, 
in agreement with recent experiments. This scenario is 
quite different from the conventional picture of a super- 
conductor, where the retardation effects are believed to 
play the central role in reducing fi* . It explains quite nat- 
urally the anomalous pressure dependence of Tc found for 
CsaCeo and the absence of a strong couphng to the al- 
kali phonons. It also predicts that the coupling to the Ag 
phonons is strongly reduced by screening effects. Finally, 
our results let us understand the surprising fact that Tc 
peaks for systems close to the Mott transition, where the 
density of states is large, but the screening has not yet 
started to become inefficient. 
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